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We investigate the dynamics of pairs of drops in microfluidic ladder networks with slanted bypasses, which 
break the fore-aft structural symmetry. Our analytical results indicate that unlike symmetric ladder net- 
works, structural asymmetry introduced by a single slanted bypass can be used to modulate the relative drop 
spacing, enabling them to contract, synchronize, expand, or even flip at the ladder exit. Our experiments 
conflrm all these behaviors predicted by theory. Numerical analysis further shows that while ladder networks 
containing several identical bypasses are limited to nearly linear transformation of input delay between drops, 
combination of forward and backward slant bypasses can cause signiflcant non-linear transformation enabling 
coding and decoding of input delays. 



Understanding the spatiotemporal dynamics of con- 
flned immiscible plugs in interconnected fluidic paths 
is essential for applications ranging from lab-on-chip 
technological'^ to physiological flows^ to porous media 
flows'^. The traflic of drops or bubbles in even simple 
networks such as bifurcating channels can be astonish- 
ingly complex due to collective hydrodynamic resistive 
interactions in the branches'^^. Although such intricate 
dynamics, in the case of lab-on-chip applications, make 
device design challenging, the collective behaviors can be 
harnessed t o pe rform useful tasks such as droplet sorting 
and storage^'^. 

Recently, the collective dynamics between pairs of 
drops have been harnessed in the so-called microfluidic 
ladder networks (MLNs) to control their relative drop 
spacing^ ^. In MLNs, two droplet-carrying parallel chan- 
nels are connected by narrow bypass channels through 
which the motion of drops is forbidden but the carrier 
fluid can leak. Current versions of ladders have fore-aft 
structural symmetry due to equally-spaced vertical by- 
passes. Such symmetric ladders are limited in function- 
ality because the distance between pairs of drops have 
been shown to decrease at the exit only for constant in- 
let flo.^™!. 

Since flexible manipulation of drop spacing 
in networks is crucial for passively regulating a variety of 
tasks including drop coalescence^^, detection and stor- 
age, there is a need to design microfluidic ladders with 
multiple functionalities. 

From a fundamental perspective, the dynamics of 
drops in MLNs is distinct compared to the widely-studied 
microfluidic loops^^^^^^. In loops, drops at junctions 
choose a given branch. These discrete choices make such 
systems non-linear, enabling coding and decoding of in- 
put signals. Since drops do not typically make decisions 
at the bypass junctions in ladders, an open question is: 
is it possible to design microlfuidic ladders that yield sig- 
niflcant non-linear transformation of input signal? 

In this Letter, we study the dynamics of spacing be- 
tween drop pairs in MLNs with slanted bypasses. We 
flnd that because the slant breaks the fore-aft symmetry, 
it provides signiflcantly more control over drop spacing 
than symmetric MLNs. We also discover that inclusion 



of slanted bypasses in ladders can non-linearly transform 
the initial delay between drops. These advanced capa- 
bilities arise because slanted bypasses flexibly manipu- 
late (i) the locations in the channels where drop velocity 
changes occur and (ii) the time drops spend with by- 
passes between them. 
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FIG. 1. (a) Ladder with a vertical bypass (b) Three dis- 
tinct configurations are possible when a pair of drops traverse 
through a symmetric ladder network; Black objects represent 
drops. Full and dashed lines denote transport and bypass 
channels respectively. Arrows show flow direction. 

As shown in Fig. [T] the key framework for quantifying 
drop spacing comes from understanding the variation in 
relative velocity {u) between drops in the top and bottom 
channel as they cross nodes in the network. Consider the 
simple case of a symmetric MLN with one vertical bypass. 
As shown in Fig. [TJb) left, when two drops driven by a 
constant flow rate enter the ladder with an initial separa- 
tion (Axin)^ they maintain the same separation, as = 0. 
When the leading drop crosses the node, a new conflgu- 
ration (see Fig. [TJb) center) is reached, and the relative 
velocity changes as fluid leaks into the bypass, i.e. u < 0, 
causing a contraction in inter-drop distance. When both 
drops cross the bypass, u = again as the pressure drop 
across the bypass is zero due to equal downstream branch 
resistances (see Fig. [TJb) right). Thus, the outlet drop 
spacing, {Ax out) ^ is less than the initial separation and is 
given by Ax out = Ax^n+uAT, where AT is the duration 
drops remain in the particular conflguration of Fig. [ijb) 
center. 

In networks with many bypasses, more drop conflgura- 
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tions and relative velocity changes are possible than that 
of Fig. [TJb) as drops traverse through multiple nodes. In 
general, we find that 

p 

Ax out = AXin + ^ Uj^Tj (1) 

j = l 

where p is the number of distinct configurations of 
droplets occurring in the network, uj and ATj are the 
associated relative velocities and time periods. 

The degree of separation achieved at the ladder exit 
depends on the contribution of the summation term in 
Eqn ([T]) both with respect to magnitude and sign. The 
strength of this contribution is modulated by the spe- 
cific architecture of the ladder network. To conipute this 
contribution, we use resistive network modelingi^, where 
we assume each drop is a point object with the same 
hydrodynamic resistance (Rd)- Since the drop velocity 
(V) is linearly dependent on liquid flow rate we 
have V = PQ/S , where S is the channel cross-sectional 
area and < /3 < 2. To preserve the relative separa- 
tion when drops leave the ladder network, we choose the 
downstream channel resistances to be equal. 

We begin by discussing the effect of a single slanted 
bypass in regulating the dynamics of drop spacing. Fig. 
[2] (a, b) shows a representative ladder network with a 
slanted bypass. In contrast to the vertical bypass, a new 
control parameter, AL is needed to describe two possible 
structural configurations — backward slant for AL < 
and forward slant for AL > 0. 

To determine the drop spacing at the exit due to the 
slanted bypass, we identify the relative velocities that are 
non-zero and the corresponding durations as prescribed 
by Eqn ([T]). Similar to the vertical bypass in Fig. [l] non- 
zero u occurs only after one of the droplets crosses a node. 
Solving the hydrodynamic circuit for this drop configura- 
tion analytically, we obtain u = j3Q/ S ■ {Rd/{Rb + 2i?e + 
Rd)) and AT = {\Axin - AL\)S/{pQ), where Rj, and Re 
are the bypass and exit channel resistance respectively. 
Thus, Eqn ([T]) for the case of an MLN with a slanted 
bypass transforms to 

AXout = AXin - M{AXin - AL) (2) 

where M = Rd/{Rb + 2i?e + Rd) and < M < 1. Note 
in Eqn Axin > corresponds to the top drop leading 
over the bottom drop in the ladder. 

Remarkably, Eqn ([2| captures several dynamical 
regimes emerging from structural asymmetry due to a 
slanted bypass as illustrated in Fig. [2|c). For the partic- 
ular case of a vertical bypass (AL = 0), Eqn Q reveals 
that drops can only undergo contraction (see Fig. [ija)). 
Perfect synchronization of drop pairs, i.e. Axout=^^ is 
difficult to achieve with a single vertical bypass, as it 
requires M to be unity. 

In contrast to the vertical bypass, we find that the 
backward slant, where AL < 0, yields flexible control 
over drop spacing as illustrated in Fig. [2|ii). For large 
input drop spacing, Axin > the pairs at the exit 
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FIG. 2. Ladder networks with (a) backward slant (b) forward 
slant; (c) Dynamical regimes due to drop traffic in MLNs with 
a single (i) vertical bypass (ii) backward slant and (iii) forward 
slant. M = 0.45. 

undergo contraction for the same reason as in the vertical 
bypass. Perfect synchronization can also be realized with 
just a single backward slant when Ax in = More- 
over, when Axin < "^^j, a new regime emerges that we 
refer to as flipping. We observe that the leading droplet 
is initially ahead of the lagging droplet. However, when 
the leading drop crosses the bypass first, its velocity is 
reduced and the lagging drop has sufficient duration to 
catch up and overtake it. Thus, the flipping behavior 
yields Axout < as shown in Fig. [2|ii). 

Similar to the backward slant, the forward slant (where 
AL > 0) provides additional means of control. Interest- 
ingly, in contrast to the backward slant, the behavioral 
transitions depend only on the value of input drop spac- 
ing relative to AL. When Ax in < AL, the lagging drop 
crosses the bypass first, resulting in expansion as high- 
lighted in Fig. |2|iii). Alternatively, if Axin > AL, the 
leading droplet crosses the bypass first, resulting in con- 
traction. If Axin = AL, then both drops cross the bypass 
simultaneously, and the input spacing is preserved. Thus 
forward slant allows drop pairs to expand, contract or 
remain unchanged. 

It would be misleading to visualize the slanted ladder 
network as being equivalent to a vertical bypass network 
with the addition of AL to the inlet spacing as we have 
already shown that vertical bypasses can only reduce the 
drop spacing for a fixed inlet flow. We assumed the exit 
channels to have equal lengths to simplify the analysis; 
however, the analysis presented in this letter extends triv- 
ially to unequal exit lengths also. 

To confirm the different behaviors predicted by our 
theory, we sought to construct an MLN with a single 
slant. However, the design space is large, requiring op- 
timization of upstream and downstream transport chan- 
nel resistance, slant resistance, slant slope, and hydrody- 
namic resistance of drops, which itself is a con iplex func- 
tion of the flow conditions and fluid propertie^^^. For- 
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FIG. 3. Experimental confirmation of tfie dynamic befiaviors 
in MLNs: (a) Snapsfiots sfiowing syncfironization of a droplet 
pair. Time interval between images is 0.08s; (b) Drop spac- 
ing as a function of position in the ladder (|AL| = 500/im). 
Continuous phase was hexadecane, and dispersed phase was 
aqueous dye solution. The transport channels are lOO/xm wide 
and tall. 



ladders with n bypasses, the maximum number of con- 
figurations where u is non-zero is n(n + l)/2. However, 
all these configurations need not be realized for a given 
Axin^ which makes the theoretical analysis complex. We 
therefore used resistive network based simulationJI^ to 
fully quantify the exit drop spacing for arbitrary input 
delay. As shown in Fig. |4ja), we find that additional ver- 
tical bypasses simply amplify the contraction effect due 
to a single bypass, i.e., slope decreases with increasing 
number of bypasses. A similar outcome also holds, as 
shown in Fig. |4]^b), for the particular case of an MLN 
with forward slants. We also find that the largest change 
in drop spacing occurs in the first few bypasses. Ladders 
with multi-bypasses could therefore be useful to maintain 
the same behavior as their single bypass counterparts, 
while dampening the effect of small fluctuations in input 
drop spacing. 
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tunately, insights from Eqn Q reduce the search space. 
According to Fig. [3) the backward slant is the best can- 
didate to achieve maximum contraction, perfect synchro- 
nization and flipping. Moreover, Eqn ([2) reveals that if 
Axin < (.i.e., the top drop is lagging behind the bot- 
tom drop in the ladder), then the forward slant becomes 
a backward slant, allowing access to the expansion regime 
as well as the condition where the input separation does 
not change. Thus, we chose the backward slant to test 
our predictions. 

We incorporated two flow-focusing drop generators in 
polydimethyl(siloxane) devices to introduce drops at a 
constant flow rate into the ladder. To amplify the effects 
produced by the backward slant, we maximized the value 
of M by minimizing Rij and Re. For example, as shown 
in Fig. [3]^ a), the bypass has an enlarged mid-section to 
minimize i?5. We also ensured and to he ^ 0{Rd) . 
In Fig. [sjb), we show experimental curves corresponding 
to each of the dynamical regimes predicted by our theory. 
We find the drop spacing to be relatively unchanged when 
the drop pair is before or after the bypass. However, as 
pairs of drops cross the bypass section, their spacing may 
expand, contract, flip, synchronize or remain unchanged 
depending on the initial separation. By comparing the 
data of Fig. |3] with the expression for we estimate 
Rd ^ 1.2kg /mm^s; M ^ 0.2. Taken together, our ana- 
lytical results and experiments suggest that MLNs with 
slanted bypasses provide greater flexibility than ladders 
with vertical bypasses. 

To understand the flexibility due to additional by- 
passes, we investigated ladders containing several iden- 
tical bypasses. In contrast to the single bypass case, in 



FIG. 4. Ladder networks with multiple identical bypasses: 
(a) Vertical bypasses (b) Forward slants. Rd/Rb = 3 
Rd/RAL 22,Rd/Re 2,Rd 1.5kg/mm^s and /3 1.4. 



Our analysis of MLN designs with identical bypasses 
has shown that at small input delay, drop spacing may ei- 
ther contract or expand, while it always contracts at large 
input delay (c.f. Fig. [4|. The contraction at large input 
delay is expected because the leading drop is the first 
to cross the bypass and slow down. Initially, it appears 
that expansion is not possible at large input delays. To 
further probe this notion, we developed an evolutionary 
algorithm^^ to search for ladder designs containing any 
combination of slant and/or vertical bypasses that might 
be capable of contraction at low input spacing and expan- 
sion at large input spacing. Our search strategy revealed 
that such networks do exist, and an example is shown in 
Fig. (sja). The dynamics of drop spacing in this network 
cannot be rationalized from mere addition of functional- 
ities of the single bypasses shown in Fig. |4] Instead, we 
find that the first five bypasses collectively cause contrac- 
tion of drop spacing, while the last two bypasses cause 
drops to expand. However, the relative magnitudes of 
contraction and expansion from these sets of bypasses 
depends on the input drop spacing. We observe that at 
small input delays, the first five bypasses dominate, re- 
sulting in contraction behavior (see Fig. [5|a)), whereas 
at large input delays, the last two bypasses dominate, 
yielding expansion. 

A striking observation from Fig. [sja) is that the curve 
is significantly nonlinear compared to the almost lin- 
ear dependence observed in ladders with identical by- 
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FIG. 5. (a) Nonlinear output delay in a ladder network 
containing a mixed combination of slanted and vertical by- 
passes (structure of the network is shown in inset) (b) En- 
coding and decoding signal using ladder network Rd/Rh — 3 

Rd/RAL = 2.6, {Re2) to {Rel) is 1 .01 = 20. Rd = 

l.^kg/mrn^s and (3 — 1.4. 



passes. This result is significant because it implies non- 
linear transformation of input delay without any droplet 
decision-making at bypass junctions. A unique conse- 
quence of this nonlinear transformation is the capability 
to encode and decode input delays as shown in Fig. [sjb) 
where the entrance delay between pairs of drops repre- 
sents the input signal and the system of bypasses repre- 
sents the encoder and decoder. First consider the curves 
in Fig. [3] where the input signal is 'scrambled' in the 
bypass section, but does not revert to its original value 
and therefore is not decoded. In striking contrast, we 
find in Fig. [sja), that input signal gets encoded and de- 
coded at two different values of Ax^^ = 4.5, 7(AL), where 
Axin = ^Xout' At these input delays, we find that the 
above-discussed contraction and expansion effects intro- 
duced by sets of bypasses negate each other. Thus, our 
results highlight a new route to code and decode signals 
compared to earlier studies that use drop decision making 
events in network^. The encoder/decoder illustration in 
this letter is just an exemplar demonstration of the po- 
tential of such networks in realizing complex behavior 
without any active components. 

Given that ladders with mixed combinations of by- 
passes can display nonlinear behavior, we ask what de- 
gree of nonlinear ity can be achieved in ladders. A qualita- 
tive indication can be obtained by considering reversibil- 
ity in ladder networks. Reversibility implies that the orig- 
inal input delay is recovered when the flow is reversed. 
Ladder networks are reversible because of the absence of 
decision-making events^^. This reversibility criterion de- 



mands that the relationship between input and output 
delays remain bijective. In addition, this functional rela- 
tionship cannot have maxima, minima, or saddle points 
because it has to be strictly monotonic. Thus, we believe 
reversibility imposes bounds on the degree of nonlinearity 
that can be achieved with microfluidic ladders. 

In summary, we observe that MLNs with mixed by- 
passes display rich dynamics in drop spacing as well 
as nonlinear behavior. Such ladder networks in fact 
also exist in natural systems including leaf venation^^, 
microvasculaturd^^ and neural systems^^ Interestingly, 
the nonlinear shape of Fig. [sja), resembles that of the 
widely observed sigmoid function, which is essential for 
coding/decoding neural signal^, as it produces an in- 
vertible map. Finally, the framework described here can 
be expanded to explore not only pairs of drops, but also 
trains of drops to further probe collective hydrodynamics 
in drop-based microfluidic networks. 
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